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The energy-based stochastic extension of the Schrodinger equation is a rather special nonlinear 
stochastic differential equation on Ifilbert space, involving a single free parameter, that has been 
shown to be very useful for modelling the phenomenon of quantum state reduction. Here we con- 
struct a general closed form solution to this equation, for any given initial condition, in terms of a 
random variable representing the terminal value of the energy and an independent Brownian motion. 
The solution is essentially algebraic in character, involving no integration, and is thus suitable as a 
basis for efficient simulation studies of state reduction in complex systems. 
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The standard energy-based stochastic extension of the 
Schrodinger equation is given by the following stochastic 
differential equation: 

dIV't) = -iH\y^t)dt - la^{H - Htf\i^t)dt 

+ ^a{H - Ht)\^t)dWt, (1) 

with initial condition \ipo). Here \ipt) is the state vector 
at time t, H is the Hamiltonian operator, Wt denotes a 
one-dimensional Brownian motion, and 



(2) 



is the expectation of H in the state |i/'t)- The parameter 
a, which has the units a ~ [energy] ~ ^ [time] 



-1/2 



gov- 



erns the characteristic timescale th associated with the 
collapse of the wave function induced by (|l|). This is 
given by Tfl = 1/a^Vo, where Vq is the initial value of the 
squared energy uncertainty, which at time t is 



Vt 



(3) 



The stochastic equation (|l|) provides perhaps the sim- 
plest known physically plausible model for state vector 
reduction in quantum mechanics [|l],H . Although its prop- 
erties have been studied extensively, it has hitherto been 
necessary to resort to numerical methods to solve (|^). 
The purpose of this article is to present an analytic so- 
lution for the dynamics of {ipt}- Apart from its use as 
a means for generating a general solution to a nonlinear 
problem in quantum state dynamics, the method we pro- 
pose also sheds new light on the nature of quantum prob- 
ability and some of the issues associated with the flow of 
information when quantum measurements are made. 

We begin with a brief overview of the stochastic frame- 
work implicit in the extended Schrodinger dynamics 
given by equation (|^). We follow closely here the analysis 
presented in a. Specifically, we introduce first the key 



notions of filtration, conditional expectation, martingale, 
and potential. We then demonstrate that the conditional 
expectation ( [l0|) gives rise to the energy expectation pro- 
cess (||) . As a consequence, we are led to simple analytic 
expressions for the energy (^5|) and the state vector ( [2^ ) 
in terms of a pair of underlying state variables. These 
results open up the possibility of efficiently simulating 
the reduction process for a variety of models. Finally 
we illustrate the practical advantages of our method by 
analysing in some detail the timescale associated with 
the reduction process in the case of a two-state system. 

The dynamics of |V't) are defined on a probability space 
(f},J^,P) with filtration JP* (0 < t < oo). Here il is the 
sample space over which J-" is a cr-field of open sets upon 
which the probability measure P is defined. 

The filtration represents the information available at 
time t. More specifically, a filtration of is a collec- 
tion {J-t} of (T-subfields of T such that s < t implies 
C J-f- Given a random variable X on {^l,T,P) 
for which E,[X] exists, we write E[A"|J^t] for the condi- 
tional expectation of X with respect to the cr-subfield 
J-t C T. Intuitively, conditioning with respect to Tt 
means giving the information available up to time t. The 
nesting J^s C J^t for s < t thus gives rise to a notion 
of causality. For convenience, we use the abbreviation 
Et[A'] = E[Ar|J^t] when the choice of {Tt} is understood. 
The conditional expectation E([Ar] satisfies: (i) the tower 
property E4Et[A:]] = Es[X] for s < t; and (u) the law of 
total probability E[Et[A:]] = E[A:]. If Et[X] ^ X we say 
that X is JFj-measuraWe. 

The conditional expectation operation allows us to in- 
troduce the concept of a martingale, the stochastic ana- 
logue of a conserved quantity. A process Xt is said to be 
an J- 1- marling ale if E[|Art|] < oo and Es[Xj] = Xg for all 
Q < s < t < oo. In other words, Xt is an J^j-martingale 
if it is integrable and if its conditional expectation, given 
information up to time s, is the value X^ of the pro- 
cess at that time. If the filtration is fixed, then we can 
simply speak of a martingale without further qualifica- 
tion. There are circumstances, however, where more than 
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one filtration can enter a problem, and then we have to 
specify with respect to which filtration the martingale 
property holds. 

For a concise mathematical representation of the state 
reduction process, we also require the concepts of su- 
permartingale and potential. A process Xt is an !Ft- 
supermartingale if E[|Xt|] < oo and Es[Xt] < Xg for all 
< s < t < oo. Intuitively, a supermartingale is on aver- 
age a nonincreasing process. A positive supermartingale 
Xt with the property E[X(] — > as f ^ oo is called a 
potential. 

The filtration Tt with respect to which stochastic ex- 
tension of the Schrodinger equation (|l|) is defined is gen- 
erated in a standard way by the Wiener process Wt- 
We signify this by writing Tt — . It is straightfor- 
ward to verify that the Hamiltonian process Ht is an 
^f^-martingale, and that the variance process Vt is an 
jFj^-supermartingale. That is to say, Es[i/t] = Hg, and 
]Es[V(] < Vs- These relations can be deduced by applying 
Ito's lemma to (|]) and (^, from which wc infer that 



Vt=M{H^ 



(7) 



Ht = Ho + G VsdW, 



and 



Vt = 



vMs 



(4) 



(5) 



Here /3t = {'(pt\{H — HtY\tpt) / {'4^t\'4't) is the skewness of 
the energy distribution at time t. The martingale and the 
supermartingale relations then follow as a consequence of 
elementary properties of the stochastic integrals appear- 
ing in d) and ^). 

In the case of the ordinary Schrodinger equation with a 
time-independent Hamiltonian, the energy process (|^) is 
constant. This is usually interpreted as the quantum me- 
chanical expression of an energy conservation law. How- 
ever, if a system is in an indefinite state of energy then 
it is not clear a priori what is meant by energy conser- 
vation. The martingale condition Es[i?t] — Hg can be 
interpreted as a generalised energy conservation law ap- 
plicable in such circumstances. In particular, it implies 
that once the state reduction has occurred, the proba- 
bilistic average of the outcome for the energy must equal 
the initial expectation. 

The supermartingale property satisfied by Vt on the 
other hand is the essence of what is meant by a reduction 
process. In fact, it follows from equation (^) that the 
asymptotic behaviour of Vt is given by limt^oo E [Vt] = 0. 
In other words, the variance process for the energy is 
a potential. Writing Hoo = Hq + cr Jg VtdWt for the 
random terminal value of the energy, one can prove |^] 
as a consequence of (^) and (||) that 



Ht=Et[H^] 



(6) 



That is to say, Ht and Vt are respectively the T}^- 
conditional mean and variance of i/oo- 

With these facts in hand, we now present a method for 
obtaining a general solution to the stochastic equation 
(|l|). The setup is as follows. We denote by Ei {i = 
1,2,...) the eigenvalues of the Hamiltonian of a given 
quantum system, and write 



l(^o|V'.)P 

(V'olV'o)(V'j|V'i)' 



(8) 



for the transition probability from the given initial state 
iV'o) to the eigenstate with energy Ei. If the spec- 
trum of H is degenerate, then denotes the Liiders 
state, i.e. the projection of |i/'o) onto the linear subspace 
of states corresponding to the eigenvalue Ei . 

Now let the probability space (il, P) be given, and 
on it specify a random variable H that takes the values 
Ei with probabilities iTi. We also assume that (r2,JF, P) 
comes equipped with a filtration Qt with respect to which 
a standard Brownian motion Bt is specified, and that H 
and Bt are independent. We assign no a priori physical 
significance to H and Bt, which are introduced as an 
ansatz for obtaining a solution for (^. 

We now define a random process S^t , which we shall call 
the signal process, according to the scheme 



^t = am + Bt 



(9) 



Intuitively, one can think of as giving a 'noisy' rep- 
resentation of the information encoded in the random 
variable H . 

We let {J-f } denote the filtration generated by the pro- 
cess ^t, i.e. the information generated by as time pro- 
gresses, and consider the conditional expectation 



Ht=E \h\T^ 



(10) 



and that 



Clearly, Tt C Qt since knowledge of H together with 
{Bs}o<s<t imphes knowledge of {Cs}o<s<t, although the 
converse is not the case. 

The significance of the JF^^-martingale Ht is that it 
represents the 'best estimate' for the value of H given 
the history of the signal process from time up to 
time t. More precisely, an ^j^-measurable random vari- 
able Yt minimises the expectation of the squared devi- 
ation of H from Yt, given T^, iff Yt = E[H\T^]. This 
follows, by a variational argument, from the relation 
- Ytf\:Ff] = E[i/2|4] „ 2YtE[H\:F^] + Y,\ 

We proceed to establish the remarkable fact that the 
process Ht defined by is statistically indistinguish- 
able from the energy process associated with the 
stochastic extension of the Schrodinger equation 

The argument goes as follows. First, because is a 
Markov process satisfying limt^oo t^^^t = H , we have 
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E 



(11) 



In other words, to determine the conditional expectation 
of H given {^s}o<s<t it suffices to condition on alone. 

To calculate E[i?|^t], we require a version of the Bayes 
formula applicable when we consider the probability of 
a discrete random variable conditioned on the value of a 
continuous random variable. In particular, 



(12) 



where tt^ = V[H — Ei). Here p{£,t\H — Ei) denotes the 
conditional probability density for the continuous ran- 
dom variable given that H = Ei. Since Bf is a stan- 
dard Brownian motion, the conditional density for is 



pi^tlH^E,) 



1 



exp ( -{^t~cjE,t)' 



(13) 



It follows from the Bayes law (|12|) that the conditional 
probability for the random variable H is 



\H^E,\it 



TTi exp {(jEiS^t - \(y'^E'lt) 



(14) 



Therefore, we deduce that the conditional expectation of 
H given S,t is 



E,¥{H = E,\^t) 



J2, ^^^E^ exp {aE.^t - ^a^Eft) 
TT, exp {<jE,^t - h<^^Eft) 



(15) 



In order to show that Ht is the energy process of the 
given quantum system, one further key result is required: 
namely, that the process Wt defined by 



(16) 



is an -Brownian motion. To verify this, it suffices, 
by virtue of Levy's characterisation of Brownian motion 
[Ql, to demonstrate (a) that Wt is an JFj^-martingale, and 
(6) that {dWty = dt. To verify (b) we note that (|) 
implies d^f — aHdt + dBt, and thus (d^t)^ ~ dt. On the 
other hand, (|l^ ) implies that dWt = d^t — crHtdt, and 
hence {dWt)^ = (dCt)^. To estabhsh (a), let (|l|) define a 
function -ff (^, t) of two variables such that Ht — H{^t, t): 



7T,E, exp {aE,^ - ^(J^Eft) 
E,^^exp {aE,^-^^a^Eft) 



(17) 



Then applying Ito's lemma and using the relation 
(dft)^ = dt, we obtain 

dHt - {dtH{^t,t) + hJl^(^t,t))dt + d^H{^t,mt (18) 



where dtH{£_t,t) denotes dH{£^,t)/dt valued at ^ = ^t, 
and so on. A short calculation making use of (0) shows 
that d^Hi^t,t) ^aV{^t,t) and dtH{^t,t) + ^dlH{^t,t) = 
-(^'^V{it,t)H{^t,t), where the function F(^,t) is 



V{i.t) 



j:,Tri{Ei-Hii,t))^e^p{aE,^- 
TT, exp (aE,^ - ^a^Eft 



f —■ (19) 



It follows, in particular, that -ff (^, t) is monotonic in f 
for any given value of t. Substituting these results into 
(p^), we infer that dHt = (^V {if,t)dWt. However, we 
know that Ht is an J-f -martingale, and therefore we con- 
clude that Wt is also an JF^^-martingale, and that estab- 
lishes (a). We thus deduce that Wt is an .F|-Brownian 
motion, as claimed. We call Wt the innovation process 
associated with Ht. The significance of Wt is that the 
process defined in (^ satisfies the diffusion equation 
d^t = crHtdt + dT4^t, where Ht = i/(Ct, t). 

Now let \4>o) be the initial normalised state vector of 
the quantum system, and let Pi denote for each value 
of i the projection operator onto the Hilbert subspace 
corresponding to the energy eigenvalue Ei. We let \ijji) = 

— 1/2 " 

TTj Pi\ipo) denote the Liiders state corresponding to Ei, 
and write Ha = P (i? = i?i|Ct) for the process defined by 
(p^. Then, we can verify that 



1^*) = E 



(20) 



satisfies the stochastic extension of the Schrodinger equa- 
tion (|l|) with the given initial condition. In particular, 
by applying Ito's lemma to ( [T^ ) and using the diffusion 
equation satisfied by ^t we obtain 



dn,t = (7{E, - Ht)UudWt. 



(21) 



With another application of Ito's lemma we deduce that 



dull' = -l<^HE^ - Ht)^Ul['dt + ^<j(E, - Ht)ull^dWt. 
A short calculation then shows that (gO|) satisfies (^, and 
that the expectation of the operator H in the state \tpt) 
is the process (psf). 

Summing up, the stochastic equation (|l|) can be solved 
as follows. We let H he a random variable taking values 
Ei with the probabilities tt, defined by (^), or equiva- 

lently tt^^ = (i/'o|A|i/'o)/(V'o|'0o)- Letting Bt denote an 
independent Brownian motion, we set ^t = (jHt + Bt. 
The solution of (|^) is then given by 

1^.^^ _ %/^exp {-iE^t + \nE,^t - l<J^Eft) \^,) ^^^^ 



(E.TT.exp (aES-^cT^Eft)) 



1/2 



where the J^^^-Brownian motion Wt driving \ipt) in (0) is 
given by (p^. In particular, by use of (p^), the expression 
(|l|) for Ht follows at once since {ipilHltpj) = EiSij. 

The fact that (|l^) is indeed a reduction process for 
the energy can be verified directly as follows. Suppose, 
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in a particular realisation of the process Ht , the random 
variable H takes the value Ej for some choice of the index 
J. Setting LOij = Ei ~ Ej and writing ~ uEjt + Bt, we 
have, for the corresponding realisation of Ht, 



Ht = 



n^E^ exp {aE.Bt - \(t^E,{E, - 2Ej)t) 
Y,i T^i exp (aEiBt - ^a^EiUJijt) 
''^j^j + '^i^i exp {auoijBt ~ ^a^ujfjt) 
'^3 + '^i exp {aujijBt - ^a^Lofjt) 



(23) 



where = However, the exponential martin- 

gale Mijt defined for i ^ j by 



Mijt = exp {aUijBt - icr^cj^^^t) 



(24) 



that appears in expression (g3|) has the property: 
limt^oo P {Mijt > 0) = 0. Hence from 



Ht = 



ijt 



'^3 + Ej 



(25) 



ijt 



we see that Ht converges to the value Ej with probability 
one. A similar argument allows us to verify that if iJ = 
Ej then for each value of i we have limt_,oo ^it — l{i=j}, 
where 1 denotes the indicator function, which shows that 
\4't) converges to the Liiders state corresponding to the 
energy eigenvalue j with probability one ||^ . 

Therefore, we see that the random variable H can be 
identified with the terminal value Hoc of the energy pro- 
cess. The fact that H is not JF^^-measurable for i < oo 
indicates that the 'true value' of H is 'hidden' until the 
reduction process is complete. On a related interpreta- 
tional point we note that in stochastic models for quan- 
tum state reduction it is sometimes assumed that the 
driving process Wt is in some way 'external' to the quan- 
tum system, representing, e.g., a noisy environmental 
coupling. This assumption, however, is unnecessary: as 
far as the flow of information is concerned, we have 



w 



— t ^ t ^ ■'t I 



(26) 



and it is thus perfectly consistent to regard the innovation 
process Wt as being endogenous. 

The advantage of the expressions (|l|) and (|2|) is that 
Ht and |'0t) are expressed algebraically in terms of the un- 
derlying random variable H and the independent Brow- 
nian motion Bt . These can be thought of as representing 
independent state variables for the reduction dynamics. 
As a consequence, we are able to investigate properties of 
the reduction process (|^) directly without having to re- 
sort to numerical integration. In particular, by use of ( |2^ ) 
a numerical simulation of the state reduction of rather 
complex quantum systems is feasible, including cases for 
which the Hamiltonian has a nondiscrete spectrum. 

In conclusion let us analyse now in detail the timescale 
associated with the reduction process. For simplicity, we 
consider a two-state system with energy levels Ei and 



E2. The initial state is given by |'0o), and the transition 
probabilities to the energy eigenstates \Ei) and \E2) are 
given by tti and ■n2- 

Suppose a measurement of the energy is made, and we 
condition on the outcome of the measurement being Ei. 
In that case, according to (p5|), we have 



Ht - 



■niEi + 7r2E2M2it 

TTl -I- TT2M21t 



(27) 



where Af2it = exp ((Ta>2iSt — ^cr^o^l^^i) . Writing /3 = 
j<j'^uj2i for the parameter that determines the character- 
istic rate of reduction, we can work out the probability 
that M2it < e~" for some value of n. Since Bt is nor- 
mally distributed with zero mean and variance t, we find 
that 



P (M2it < e-") =P(^Bt< - \r^'^nj 

= N (iptf^ - in(/3t)-i/2) (28) 

where N{x) is the standard normal distribution function. 

Therefore, for example, we see that provided t > bra, 
we have P (M2it < e~^°) > i, where tr = 1/(3. In par- 
ticular, as Ht draws near Ei we have the relation 



Ht-Ei 



^{E2 



Ei)M2if. 



(29) 



Thus, after only a relatively few multiples of the char- 
acteristic reduction timescale, the amount by which Ht 
differs from Ei will typically be reduced to a tiny fraction 
of the energy difference E2 — Ei. 
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